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ABSTRACT 

We present photometric and spectroscopic data on the early-type binary MP Centauri. 
The photometric data are analyzed simultaneously with radial velocities to derive 
preliminary absolute dimensions for the binary components. Analysis of the spectra 
shows that the stars rotate synchronously and that the line of sight to the system 
crosses two kinematically sharp and well-separated interstellar reddening sources. It 
is shown that MP Cen consists of a B3 primary with Mi = 11.4 ± 0.4Mq, i?i = 7.7 ± 
0.1i?Q and a lobe-filling B6-B7 secondary with M 2 = 4.4±0.2Mq, R 2 = 6.6±0.1i?Q 

Key words: binaries: eclipsing - binaries: spectroscopic - stars: individual (MP Cen) 
- stars: evolved. 


1 INTRODUCTION 

MP Centauri (HD 308976) is an eclipsing binary, listed in 
the HD catalog as having a B3 spectral type. Its orbital pe¬ 
riod of 2.99 days makes it a very difficult target to obtain 
complete light curves from a single location in one observing 
season. It is a member a rare group: massive eclipsing bina¬ 
ries with short orbital periods. It lies in the galactic plane 
(6 = 0°.08, I = 295°.01) and its eclipsing nature allows for 
the determination of its absolute parameters and distance, 
as well as the reddening along the line of sight. We have ob¬ 
tained uvby photometry over several years and the system 
was als o observed in V a nd Iq by the All Sky Automated 
Survey iPoimanskill200 j : hereafter AS AS). 

As part of a program on early-type overcontact and 
near-contact binaries, we observed MP Cen with the ESO 
La Silla 2.2m telescope and the FEROS spectrograph at high 
resolution (R=48,000) in February, 2003. Combining radial 
velocities from these spectra with the uvby photometry and 
V and Ic data from ASAS, we present the first comprehen¬ 
sive study of this evolved binary. 


* E-mail: terrell@boulder.swri.edu (DT); munari@pd.astro.it 
(UM); tomaz.zwitter@fmf.uni-lj.si (TZ); georgewolf@smsu.edu 
(GW) 


2 OBSERVATIONS 

2.1 Photometry 

GW observed MP Cen in 1982 and 1989 using the 24- 
inch (0.6 meter) telescope and single channel photometers 
equipped with uvby filters at the Mt. John Observatory at 
Lake Tekapo, New Zealand. The comparison star used was 
HD 102139, for which the SIMBAD database gives a B4 
HI spectral type. A 17 arcsec aperture was used and offset 
slightly to avoid contamination from nearby companions. 

The ASAS project has also obtained V and Ic CCD 
photometry of MP Cen. Because the goal is to survey large 
areas of the sky, the ASAS setup is not optimal for pho¬ 
tometry of stars in crowded fields like MP Cen. However, 
MP Cen is the brightest object in the measuring aperture 
and a proper treatment of the third light {I 3 ) contamination 
makes it possible to use the ASAS data. The scatter in all 
passbands was about 2%. 

2.2 Spectroscopy 

Seven ESO 2.2m -I- FEROS spectra have been secured, in 
two groups of three spectra each on 2003 Feb 23 and 24 near 
quadratures and a single spectrum on 2004 Feb 25 in pri¬ 
mary eclipse (Table 1). The wavelength range of the spectra 
extends from 3900 A to 9200 A with R = 48, 000 and a 1200 
second exposure time. The S/N is evaluated around 5850 A 
and the orbital phase is computed with the ephemeris given 
in section 13 

Our estimate of the spectral type of the primary, based 
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Table 1. ESO 2.2m + FEROS Spectra of MP Cen. 


Spectrum Number 

Date 

UT of mid-exposure 

S/N 

Orbital Phase 

1440, 1441, 1442 

2003 Feb 23 

05:40 

75 

0.28 

1552, 1553, 1554 

2003 Feb 24 

09:15 

77 

0.66 

1632 

2003 Feb 25 

06:48 

85 

0.96 


on the relative intensities of Mgll 4481, Hel 4471, Cll 4267, 
He I 4388, Hei 4009 and Hei 4026, is B3 with an estimated 
error of one spectral subclass. Because of the relative faint¬ 
ness of the secondary and its lines, it is more difficult to 
estimate its spectral type but a value in the range B6 to 
B7 seems reasonable, based on the ratio of Hei 4471 to 
Mg II 4481, and is consistent with our light curve solution. 


2.2.1 Rotational velocity 

To derive the rotational velocities we used the relation 

Viot = 42.42 X HIW — 35km sec~^ (1) 

calibrated bv iMunari fc Tomasellal ^1999^ on the width at 
half maximum (HIW) of He l 5876 A in high resolution spec¬ 
tra of O and B stars. The deconvolution of the Hei 5876 A 
line profile in the MP Cen quadrature spectra gives a HIW 
of 4.12 A for the primary, and 3.18 A for the secondary, 
which corresponds to 140 and 100 km sec“^ respectively. 

The Munari and Tomasella relation was calibrated using 
single, and thus axially symmetric, stars. In previous w ork 
on the overcontact binary TU Mus llTerrell. et al.ll^03^ . we 
found that it was necessary to correct the HIW for the dis¬ 
torted shapes of the two stars. The MP Cen primary is rea¬ 
sonably spherical so no correction is necessary, but the sec¬ 
ondary is highly distorted s o, following the procedure out¬ 
lined in iTerrell. et al.l JioO^, we find that the corrected ro¬ 
tational velocity of the MP Cen secondary is 95 km sec“^. 
The profiles, however, are not indicative of purely rotational 
broadening, with atmospheric effects and circumstellar ac¬ 
tivity probably distorting the profile. In the absence of more 
sophisticated modeling, we see no evidence of asynchronous 
rotation and thus assume synchronism in the rest of our 
analysis. 


2.2.2 Reddening 

iMunari fc Zwitt^ lll997^ demonstrated a method of using 
the interstellar Nal and Kl lines to estimate the redden¬ 
ing affecting stars, independent of the knowledge of the in¬ 
trinsic colors. The prohle of Nal D1 and D2 lines in our 
FEROS spectra of MP Cen are shown in Figure Q and are 
clearly composed of at least two components. The 7 individ¬ 
ual spectra give very consistent results and indicate a value 
of E{B — V) = 0.30. The two-Gaussian mean fit is overplot¬ 
ted on the D2 line in Figure 0 to show the accuracy of the 
ht. 

The heliocentric radial velocities of the two compo¬ 
nents are 6.23 ± 0.13 km sec”*^ and 23.57 ± 0.08 km sec“^. 
The accuracy of these velocities is not hctitious because 
nearby telluric absorptions show the same wavelength sta¬ 
bility from one spectrum to the other, and it matches the 
expected high FEROS spectrograph stability, suitable for 



Figure 1. The interstellar lines of Nal, averaged from seven spec¬ 
tra, with the two-component Gaussian fit to the D2 line used in 
the reddening determination. The similar shape of the D1 line 
shows that the shape of D2 is not an artifact of noise or atmo¬ 
spheric lines. 


extra-solar planet searches. The line of sight to MP Cen 
therefore crosses two kinematically very sharp and well sep¬ 
arated sources of reddening of purely interstellar origin since 
neither of the components shows a radial velocity change 
with orbital phase nor shares the systemic velocity of the 
binary. 

T he Tych o Bt — Vt = 0.13 ± 0.04 corresponds, accord¬ 
ing to iBesselJ ll200Cll . to a Johnson {B — V) = 0.11. We 
denote an intrinsic (i.e., unreddened) color by {B — V)' and 
an observed color hy {B — V) with subscripts 1, 2 and C 
denoting the primary component, the secondary component 
and the composite value for the binary respectively. With 
an assumed intrinsic color for the primary, {B — V)i, and 
the magnitude differences between the two components in 
each filter, Atub and Amy, from our light curve solution, 
we can compute the intrinsic composite color of the binary, 
(B-VYc, as 


{B - VYc = {B- VYi - 2.5 log 


I _|_ ]^Q(-0.4xAmfl) \ 

T^PToPF4xA77(7y ) 


Assuming a v alue of {B — VY = —0.22 for the B3 primary 
jFlowerl ll996tl and the Ams and Amv values of 1.16 and 
1.12 respectively from the light curve solution, we find {B — 
V)c = —0.21. This result, combined with the observed (B — 
V)c, yields a reddening of E{B — V) = 0.32, in excellent 
agreement with the results from analysis of the interstellar 
Nal and Kl lines. 
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heliocentric wavelength (A) 

Figure 2. The Hel 5876 profile of MP Cen for the averaged Feb 
24 spectra. The overplotted curve is the fit of two gaussians at 
the radial velocities of the two components given in Table 2. 


Table 2. Heliocentric radial velocities of MP Cen. 


TU Mus llTerrell. et 7111200.'^) . obtained in the same observ¬ 
ing run as the MP Ce n spectra, were in excell e nt agr eement 
with earlier data from lAndersen fc GrOnbechl |I197 !tI) . 

We also performed tests of our solution results by gen¬ 
erating 1000 sets of synthetic radial velocities with Gaussian 
errors based on our four velocities and their estimated er¬ 
rors. We solved each set of four velocities and found that the 
average and standard deviation of each parameter was con¬ 
sistent with the results of our full solution. In particular, the 
average value of the semi-major axis of the relative orbit, of 
critical importance in determining the absolute dimensions 
of the binary, was 21.9±0.27?q, matching exactly the result 
from the full solution- _ _ 

WD now uses iKurucd (Il99.'^l stellar atmospheres to 
model the radiation from the stars. WD can also use ei¬ 
ther binary orbital phase o r time as the independent vari¬ 
able JWilson fc Terrelil99Sll . In order to investigate possible 
changes in the orbital period, we used time as the indepen¬ 
dent variable and adjusted the orbital period (P), its first 
time derivative (P), and the reference epoch [HJDq). We 
began our light curve fitting experi ments in WD’s mode 2 , 
appropriate for detached binaries iLeung fc Wilso^ll977ll . 
but found that a semidetached configuration with the sec¬ 
ondary Hlling the lobe was required to fit the observations 


Spectra 

Phase 

Primary (km sec 

Secondary (km sec ^ 

1440-1442 

0.28 

- 97.4 ± 2.7 

+256.0 + 4.8 ^ 

1552-1554 

0.66 

+88.6 ± 2.3 

-216.9 + 3.3 j 


2.2.3 Radial velocities 

The radial velocities of the MP Cen components have been 
measured on the two groups of spectra in Table 1 obtained 
close to quadrature phases. The spectrum taken at orbital 
phase 0.96 is very close to conjunction and thus has insuf- 
hcient velocity separation between the two components to 
obtain reliable radial velocities. Line splitting at quadrature 
is wide and allows an easy and firm determination of the ra¬ 
dial velocity, as shown for Hel 5876 in Figure |21 The radial 
velocities have been obtained by Gaussian de-convolution 
of the profiles into two components, which provide a good 
overall ht to the observed line profile. The measured lines 
are Hei 4016, 4388, 4471, 5876, 6678, 7065 and half-weight 
contributions from Hy, H/3 and Ha. The means of the ra¬ 
dial velocities (and associated standard errors) derived from 
these lines are given in Table 2. 


3 DATA ANALYSIS 

We performed a simultaneous analysis of the ASAS photom¬ 
etry, our uvby photometry and our radi al velocities with the 
2003 version of th e Wilson-Devi nney iWilson fc DevinnevI 
Il97li: IWilsonlUoT^ I WilsonI 1199(1: hereafter, WD) program. 
As we had limited telescope time available to us for spec¬ 
troscopy, we could not obtain as many radial velocities as we 
would have liked and, therefore, we present our solution as 
preliminary until more extensive spectroscopic observations 
can be made. However, our spectra are of high quality and 
resolution, giving us conhdence that our results are reliable 
and certainly better than a solution based solely on the pho¬ 
tometry. We have shown previously that our velocities for 


Other parameters adjusted in the simultaneous solu- 
i were the semi-major axis of the relative orbit (a), the 
binary centre of mass radial velocity {VP, orbital inclina¬ 
tion (i), secondary mean effective temperature (T 2 ), modi- 
hed surface potential of the primary (fli), mass ratio (q), and 
the bandpass-specific luminosity of the primary (Li). Cer¬ 
tain parameters, such as the bolometric albedos and gravity 
brightening exponents, were held fixed at their expected the¬ 
oretical values. The lo garithmic limb darke ning law was used 
with coefficients from I Van Hamrn3 The mean effec¬ 

tive temperature of the primary was set to 18,750 K based 
on the B3 spectral type (iBessell. Castelli fc Ple2^1998^ . Data 
set weights were determined by the scatter of the observa¬ 
tions. 

Given the large apertures used in the ASAS photome¬ 
try and the moderately crowded held around MP Cen, we 
adjusted third light and found statistically signihcant values 
for the ASAS data but not our uvby photometry. Images of 
the MP Cen held show two companions to the northeast 
about 13” away. These stars ht within the ASAS aperture 
but were avoided in our uvby photometry. 

Table 13 shows the results of the simultaneous solution 
and Figures and |1] show the hts to the Stromgren and 
ASAS data respectively. Figure |21 shows the hts to the radial 
velocities. The light curves show some disturbances, particu¬ 
larly at the maximum preceding primary eclipse, as expected 
if a mass transfer stream from the secondary is impacting 
the primary. 

The ASAS data show that the system has V — 9.95 and 
with the third light correction from our light curve solution, 
V — 9.99 for the binary at maximum brightness. Correcting 
this value for the reddening (assuming a ratio of total to 
selective absorption of Rv ~ 3.1 and our Eb-v = 0.30), we 
hnd that the apparent V magnitude of the binary is 9.08. 
Using the V luminosity ratio L2/Li=0.36 from our light 
curve solution, this transforms into Johnson V magnitudes 
of Vi=9.41 and V 2 = 10.53. With our log luminosity of 3.8Z/q 
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Table 3. Parameters of MP Cen 



Figure 3. The Stromgren photometry and computed light curves. 



Parameter 

Value 

a 

21.9 ±0.2^0 

w 

2.7 ± 1.7 km sec-l 

i 

82°.2±0°.2 

Ti 

18, 750 K 

T 2 

12,390 ± 50 K 

Qi 

3.29 ±0.03 

4 

0.390 ± 0.008 

HJDo 

2447627.9453 ± 0.0006 

P 

2.993456 ± 0.000001 days 

P 

8.7 X lO-'^ ±6.1 X 10-1° 

L^i/ + L^2)u 

0.83 ±0.01 

+ L2)v 

0.76 ±0.01 

L\I(L\ Lpih 

0.75 ±0.01 

^\l (-^1 + -^ 2 ) 4 / 

0.74 ±0.01 

Li/(Li + L2)v 

0.74 ±0.01 

Ti/(Li -1- 7/2)/ 

0.72 ±0.01 

hiV) 

0.04 ±0.02 

h{l) 

0.05 ±0.02 

lz(uvby) 

0.00 ± 0.02 

Ri 

7.7 ±0.1 Rq 

R 2 

6.6 ± 0.1 Rq 

Ml 

11.4 ±0.4 Mq 

M 2 

4.4 ±0.2 Mq 

log LiILq 

3.8 ±0.2 

log L 2 ILQ 

3.0 ±0.2 

log gi 

3.7 ±0.2 

log 92 

3.4 ±0.2 


Quoted errors are formal l-a errors from the solution. Z 3 values 
are in units of total system light at phase 0.25. Luminosity errors 
are estimates based on an uncertainty of 2000 K in the effective 
temperature of the primary. The error contributions to the lumi¬ 
nosities due to the errors in the radii are more than an order of 
magnitude smaller. The log g values are in CGS units. 


Figure 4. The ASAS V and Ic photometry and computed light 
curves. 



Figure 5. The observed radial velocities and computed curves. 


for the primary, the distance to the primary amo unts to 
3.0 kp c nsing a bolometric correction of -1.75 from iFlowerl 
. The same calcnlation for the secondary results in the 
same distance. We therefore adopt 3.0 kpc as the distance 
to MP Cen, placing it on the far side of the Sagittarius arm. 
MP Cen lies in the galactic plane a nd we can com- 

S ur estimate of the reddening with the lNeckel fc Klarel 
extinction maps. MP Cen lies at the boundary of 
four cells (294/-1, 296/-1, 297/0 and 292/2). Given the low 
resolution of the cell boundaries, we cannot unambiguousily 
assign MP Cen to a particular cell. Our value of Eb-v=0.30 
leads to Av ~ 0.93. Looking at the individual cells, we find 

• 294/-1 distance greater than 1 kpc 
• 296/-1 the distance could be as high as 4 kpc 
• 297/0 distances from 1 to 4 kpc supported 
• 292/2 distances from 0.5 to 4 kpc supported 

Thus, the lNeckel fc Klar^ JlOSOll maps do not provide a con¬ 
clusive comparison with our derived distance, but their Fig¬ 
ure 9 supports the idea that the line of sight to MP Cen 
goes through clouds A and G and just misses clouds L, M 
and R. 
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4 CONCLUSIONS 

We have presented preliminary values for the parameters 
of MP Cen based on photometry and limited spectro¬ 
scopic data. The masses and radii of the stars, along with 
the semidetached configuration, indicate that the system is 
evolved. The radii, Ri = 7.7Rq and R 2 = 6.67?q, are con¬ 
sistent with a giant or sub-giant classification for both stars 
for the derived masses of Mi = 11.4 Mq and M 2 — 4.4Mq. 
The P we have measured implies a conservative mass trans¬ 
fer rate of 1.1 x 10“^ Mq yr~^, making it a fairly active 
system. Light curve variations from season to season, most 
likely arising from this mass transfer activity, are apparent 
in the photometric data. We also see emission features in 
our spectra but the limited phase coverage of the spectra 
makes any discussion of the circumstellar environment little 
more than speculation at this point. More extensive spec¬ 
troscopy and polarimetry will be required to fully map the 
circumstellar environment of the system. 
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